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ABSTRACT 

The long-term dynamics of Oort cloud comets are studied under the influence of both 
the radial and the vertical components of the Galactic tidal field. Sporadic dynamical 
perturbation processes are ignored, such as passing stars, since we aim to study the 
influence of just the axisymmetric Galactic tidal field on the cometary motion and 
how it changes in time. We use a model of the Galaxy with a disc, bulge and dark 
halo, and a local disc density, and disc scale length constrained to fit the best available 
observational constraints. By integrating a few million of cometary orbits over 1 Gyr, 
we calculate the time variable flux of Oort cloud comets that enter the inner Solar 
System, for the cases of a constant Galactic tidal field, and a realistically varying 
tidal field which is a function of the Sun's orbit. The applied method calculates the 
evolution of the comets by using first-order averaged mean elements. We find that the 
periodicity in the cometary flux is complicated and quasi-periodic. The amplitude of 
the variations in the flux are of order 30 %. The radial motion of the Sun is the chief 
cause of this behaviour, and should be taken into account when the Galactic influence 
on the Oort cloud comets is studied. 
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1 INTRODUCTION 

Strong observational evidence supports the idea that the in- 
ner Solar system is subject to a stead y flux of 'n ew' comets 
which originate from the 'Oort cloud' ( Oortl 19 501. The semi- 
major axes of these comets are thought to evolve under 
the influence of external forces such as the Galactic tidal 
field and passing stars. The comets evolve from a typical 
a or i g ~ 3 x 10 4 AU to either hyperbolic orbits or to larger 
binding energies, depending on the orbital evolution dur- 
ing the cometary encounters with the Solar System planets. 
During one orbital evolution, comets typically experience 
perturbations that change the comet's orbit into a short- 
period or hyperbolic orbit, which leads t o a subsequent ejec- 
tion into interstellar space l|Nurmill200lh . 

The Galactic potential exerts a force on Oort cloud 
comets, and is important f or th e stead y state flux of comets 
with a > 20000 AU, e.g. lBv3 ll983T) ; iHeisler fc Tremaind 
l|l98rjhlMatese & WhitmarI]l989F The Galactic tidal force 
has a dominant component that is perpendicular to the 
Galactic plane; the radial component is 10 times weaker 
than the perpendicular component Hcislcr fc Tremainel 
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For this reason, many studie s (e.g. iMatese et al.l 
1 19951 ; IWickramasinghe fc Napierl 12008^ have assumed that 
the radial Galactic tidal field component is negligible. 
Other studies have included the radial component, for 
m ore accurate modell i ng of cometary mot ion, such as 
in IMatese fc Whitmird l | 1996t h iBrasserl (|200lh . The radial 
component of the tide has been found to have an ef- 
fect on the long-term evolution of the comets' perihe- 
lia, on the distribution of the longitudes of the perihelion 
|Matese fc Whitmird 19961). and the origin of chaos in the 
cometary m otion llBreiter et al.ll2008l ). In recent large scale 
simulations, iRickman et al. (|2008l ~ found that a fundamen- 
tal role is played by perturbations due to passing stars on 
comets, contrary to the in vestigations during t he pr evious 
two decades, starting with IHeisler fc Tremaind (|l986). The 
stellar perturbations, do of course act together with the 
Galactic tide. 

The topic of this paper is the effect on comets due to 
the Galactic tide alone. We use a realistic model of the lo- 
cal Galaxy, which is well constrained by observations, which 
contains a disc, bulge and halo. We follow the orbit of the 
Sun in this model using recent constraints on the Solar mo- 
tion. This motion allows us to compute the change in the ver- 
tical and radial components of the Galactic tide with time, 
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and the change in the tidal force due to the radial motion 
of the Sun is fully accounted for. 

Qualitatively, the tidal effect on the cometary orbits can 
be evaluated by studying the change in angular momentum 
averaged over one orbit. The Galactic tidal force periodically 
changes the angular m omentum (J 
the Oort cloud comets ((Fernande z fc Id 



/GM Q a(l - e 2 )) of 
199 ill. The angular 



mom entum of the comets changes as (Hei sler fc Tremainel 
1 19861 ): 



dJ_ 
~dt 



5-Kpo 



L 2 (L 2 - J 2 )[l - (J 2 /J 2 )} sin2cj s 



(1) 



Here, u) g is the Galactic argument of perihelion, J z is the 
z-component of angular momentum, and is perpendicular to 
the Galactic plane, L = yfjla, po is the local mass density, 
and G is the gravitational constant. Note that this equation 
is valid for first order mean elements under the action of axi- 
ally symmetric disc tides. The periodically changing angular 
momentum causes variations mainly in perihelion distance q, 
for comets in near-parabolic orbits, since q ~ J 2 /(2GMq). 
The typical assumptions in the cometary flux calculations 
due to the Galactic tidal force suppose that the local tidal 
field is axisymmetric, perpendicular to the mid-p lane, and 
adiabatically changing (i Matese fc Whitmanl! 19921 ). 

The aim of many of these studies is to correlate the mo- 
tion of the Sun in the Galaxy with phenomena on the Earth, 
such as mass extinctions of species, the cratering record and 
climate change. An extensive r eview of this topic has been 
made recently by iBailer- Jones! l|2009h . 

For example, early studies have shown that the ages of 
well dated impact craters on Earth are n ot distributed ran- 
doml y, but that there is a possible 28 Myr l|Alvarez fc Mullet] 
Il984h or 30 ± 1 Myr periodicity in cra ter ages over the 
past 250 Myr (|Rampino fc Stothers]|l984r i. Since then, sev- 
eral authors have claimed that there is a significant peri- 
odic signal present, but the periods differ quite a lot from 
study to study. The signal is the most prominent for 40 
large , well-dated craters that are up to 250 Myr old (|Napiej 
l2006h . but the period i s difficult to m easure, with estimates 
of between 24- 26 Myr (lNapierll2q06ri. 3 Myr ilNapierllgOOrj ; 



2006), 38 



IStothersI [20061) 36 Myr dNapierl 120061; IStothersj I 
Myr llYabushital|2004j ; IWickramasinghe fc Napiei]|2008l ') and 
42 Myr (|Napie?(2006f ). 

In some studies, the reliability of the signal is ques- 
tioned altogether, based on the inaccuracy of the age esti- 
mates of the impact craters, possible biases caused by round- 
ing the ages of craters, and the small num ber of craters 
( Gri eve fc Pesonenlll996l ; iJetsu fc Pehll2000l ). 

By critically reviewing many studies that have tried to 
connect the Solar motion and periodicity in terrestrial phe- 
nomena such_as_Wodiversity, impact cratering and climate 
change, IBailer- Jones! (2009) has concluded that there is little 
evidence to support these connections. By studying the arti- 
ficial cratering data lLvvtinen et al.l (|2009T ) came to the same 
conclusion, that the reliable detection of any periodicity is 
currently impossible with the existing cratering data. 

In this study, we statistically analyse how the Galactic 
tidal force changes cometary orbits over 1 Gyr, using nu- 
merical simulations. The 1 Gyr time-scale is long enough 
to observe changes in the Galactic tide due to both the ra- 
dial and vertical motion of the Sun. A simple axisymmetric 
Galactic potential is adopted. To our knowledge, there has 



been no study to date, in which the effect of both radial 
and vertical components, in a time varying Galactic poten- 
tial (via the variation in mass density p), has been analysed 
in detail. Our purpose is to study the statistical effects of 
the complete Galactic potential to the Oort cloud comets 
in detail, especially concentrating on the comets that enter 
the Solar System (q < 30 AU). In particular, we analyse the 
differences in cometary motion for when the tidal field is 
constant, and when it varies as the Sun moves in a realistic 
orbit in a fairly realistic Galactic potential. We find that it 
is important to include the radial motion of the Sun in the 
calculations, since the local density varies significantly as 
the Sun moves towards and away from the Galactic centre. 



2 METHODS 

The method of simulation requires two, traditionally sep- 
arate, components. The first is to simulate the motion of 
the Sun around the Galaxy. The second involves the evolu- 
tion of the orbits of comets in the Oort cloud. We integrate 
the motion of the Sun in an axisymmetric Galactic poten- 
tial. The method of integ ration of the comets is described in 
iMikkola fc Nurmil ((2006). The method calculates the evolu- 
tion of the comets by using first-order averaged mean ele- 
ments. We do not include the random perturbations caused 
by the planets, since the aim is to identify the tidal effects 
of the Galactic potential on the flux of comets reaching the 
inner Solar System. 



2.1 The Galactic potential 

The Galactic potential consists of a disc, bulge, and 
dark halo, and is partially described in iGardner fc FlvnrJ 
j|2010h. The mod e l use d here differs most notably from 
IGardner fc Flvnnl l|2010l ) in the treatment of the disc, as 
well as a slightly modified dark halo. W e noticed that the 
vertic al density profile of the disc, from IGardner fc Flvnnl 
(2010), does not accurately reprod uce the observational pro- 
file from lHolmberg fc Fly nn (2004). We proceeded to modify 
our disc-model by changing a few of the model's parameters, 
and adding three more Miyamoto-Nagai potentials to the 
model, to emulate a very thin layer of gas in the disc. The 
equation for the full potential is: 



z 

GM Cl GMc 2 



$c = - - 



$1 



1\, 



E 



E 



-GM di 



-GM 9n 

^R 2 + [a dn + ^f(z 2 + b 2 )Y 



(2) 
(3) 

(4) 
(5) 

(6) 



where G is the gravitational constant, 7? is the Galacto- 
centric radius, and z is the height. The modified parameters 
are in Table [1] The vertical density profile of the model can 
be seen in Fig. [2] 
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Table 1. Parameters of the full potential. 



Property 


value 


Unit 


v h 


220 


km _1 
m s 


ni 


10 




I, 


0.45 




1, 


0.12 


kpc 


r c ± 


2.7 


kpc 


r c 2 


0.42 


kpc 


«di 


5.81 


kpc 


a d 2 


17.43 


kpc 


a d 3 


34.86 


kpc 


M Cl 


3 


10 9 Mq 


M C2 


16 


10 S Mq 


M dl 


66.06 x10 s 


Mq 


M d2 


-59.05 xlO 9 


Mq 


M d3 


22.97xl0 9 


Mq 


M gi 


18.63 x10 s 


Mq 


M g2 


-16.66x10 s 


Mq 


Mg 3 


6.48 x10 s 


Mq 



The mass density at the Sun is 0.11 Mq/pc 3 , consis- 
tent with observat io nal co nstraints (0.10 ± 0.01 Mg/pc 3 
Holmberg fc Flynnl d2000T I and 0.105 ± 0.005 M /pc 3 



Korchagin et al l l^OOSD ^ This corresponds to a nominal 



T z =83 x 10 year period for small amplitude simple har- 
monic motion in the vertical direction. The surface den- 
sity of disc matter in the model is 54.9 Mq/pc 2 , com- 
pared with a measu r ed dis c surface density of 56 ± 6M Q /pc 2 
Holm berg fc Flvrml (2004). The adopted current position of 
the Sun in the model, (R, z)q, is (8,0) kpc. The local circular 
velocity of the model is 221 km s _1 . 

2.1.1 Density and rotation curve 

Fig. Q] shows the surface density of the disc with radius (R), 
the change in density with height (z) at the Sun (7? = 8kpc) 
(Fig. [2]), and the rotation curve (Fig. [3]), as these are the 
two factors that have the most impact on the orbits of the 
comets. The disc has a scale-length of 3 kpc, and a local 
scale-hei ght of 0.24 kpc, consistent with recent measure- 
ments by Juric et all l)2008l ) using the Sloan Digital Sky Sur- 
vey (SDSS). 

The potential is axisymmetric, so we do not model ef- 
fects such as molecular clouds, spiral arms, and bubbles in 
the interstellar matter and passing stars. Our interest here 
is on the global Galactic effects of the Solar motion on the 
comets. 
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Figure 1. The surface density of the disc component as a func- 
tion of Galactocentric radius. The dashed line corresponds to an 
exponential density falloff of 3 kpc, which is a good fit to the 
model over a wide range of radii. Note that the density truncates 
strongly at 18 kpc. 
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Figure 2. Vertical density of the model, at the Sun (R = 8 kpc). 
The dotted line represents the baryonic contribution in the model, 
the dashed line the dark matter contribution, and the dashcd- 
dotted line the total density of the model at a certain height (z). 
The solid line corresponds to an exponential fit to the baryonic 
component of the model of 0.25 kpc. 



2.1.2 Motion of the Sun 

The orbit of the Sun, using the Solar motion measured by 
ISchonrichet"ai1l|2010l ). where (U, V, W)=(U.l, 12.24, 7.25) 
km s~\ and assuming (R, z)q to be (8,0) kpc, is shown 
in the (R, z)-plane in Fig. [3J U is the velocity towards the 
Galactic centre, V is the velocity along rotational direc- 
tion, and W is the velocity perpendicular to the Galactic 
plane. Properties of this orbit, such as eccentricity, radial 
and vertical period and maximum z height, are shown in 
Table [2] The eccentricity of the Solar orbit (e) is defined as: 

max -R mm 

)/{R 

max + -Rmin). 



Table 2. Eccentricity, e, maximum vertical height, z m ax, radial 
oscillation period, Tr, and vertical oscillation period, T z , of the 
S un in the adopte d pote ntial. The adopted solar motion is that 
of lSchonrich et al.l ll2010i) . 



Property 


value 


Unit 


e 


0.059 ± 0.003 




Zmax 


0.102 ± 0.006 


kpc 


Tr 


149 ± 1 


Myr 




85 ± 4 


Myr 



(R 
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Figure 3. Rotation curve for the model Milky Way (dotted), 
and the different contributions of the disc (dashed) , bulge (long 
dashed), and dark halo (solid). 




radially in the potential. Secondly, we assume that there is 
a 'constant' tidal field (i.e. the tidal field does not change 
as the Sun moves around the Galaxy). In the constant case 
the Sun moves on a flat, circular orbit. 

We integrate the Solar orbit in the 'dynamic' case for 1 
Gyr, sampling the values of the G-parameters every 100 kyr: 
these are shown in Fig. [5] The changes in G3 are dominated 
by the changes in local density during the orbit, since the 
changes in A and B over the orbit are not very large. Fig. [5] 
(top and centre panel) shows how the changing values of the 
Oort constants, A and B, affect the values of Gi and G2. 

Fig. E] shows the combined effects of the radial (R) and 
vertical motion (2) of the orbit, on G3. It is clear that G3 
increases in two situations: when the radial position is the 
closest to the Galactic centre, and when the vertical motion 
crosses the mid-plane. Due to the slightly eccentric motion 
of the orbit, it is the radial component of the motion which 
dominates the changes in local density (p(R, z)), rather than 
the vertical mot ion. This is partly bec ause we adopt the 
Solar motions of ISchonrich et all (|2010l ). which produces a 
mildly eccentric orbit for the Sun (e = 0.059 ± 0.003): it 
oscillates between Galactocentric radii of 7.9 and 8.9 kpc. 
As such, the evolution of G3 depends on the radial motion 
and the vertical motion, both being of equal magnitude (Fig. 

In the second case studied, the Sun is set on a perfectly 
circular and flat orbit, so that the local mass density does 
not change with time. For the 'constant' case, the values of 
d are the mean values from the 'dynamical' case, and have 
the values Gi = -7.00897 x 10" 16 , G 2 = 7.27613 x 10~ 16 , 
and G 3 = 4.55749 x 10~ 15 yr -2 . 



0.15 I 1 ' 1 ' ' 1 
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R(kpc) 

Figure 4. Motion of the Sun in the model potential for 1 Gyr, 
in radial R and vertical z components. The initial position of the 
Sun is (R, z)o = (8,0) kpc. 

2.2 Variation of the tidal parameters caused by 
the motion of the Sun 

Our integration method for cometary motion is based on 
the computation of secular motion of the comet in quadratic 
perturbation by the Galactic potentia l. The method used is 
presented in lMikkola fc Nurmil l|2006l ). As the comet orbits 
the Sun under the Galactic potential it experiences a force 
F per unit mass: 

F = -^^-r-G lX x-G2yy-G 3 zz, (7) 

where Gi = —(A - B)(3A + B), G 2 = (A- B) 2 , and 
G 3 = tnGp(R, z) — 2(B 2 — A 2 ). Here r is the Sun-Comet dis- 
tance, p(R, z) is the local mass density, an d G is the gravita- 
tional constant l|Heisler fc Tremaind[l986l ). A and B are the 
Oort constants, and are obtained from the Galactic model. 
Usually the radial components (x,y) are neglected, so that 
Gi = G 2 = 0. 

We will examine the effect on comets in two particular 
cases. Firstly we study what we call the 'dynamic' case, in 
which the Galactic tide changes realistically along the Solar 
orbit, for the case that the Sun oscillates both vertically and 



3 SIMULATIONS OF THE COMETS 

Observations support the idea that the in situ flux of 
new comets is roughly linear with respect to heliocentri c 
distance, up to the distance of Jupiter ( jHughes! 1200 ll ). 
Comets with highly eccentric orbits with randomis ed di- 
rectio ns of motion have a g-distribution that is flat (|0pikl 
1966). If a long period comet has a perihelion distance 
between ~ 10—15 AU, it is quickly ejected to interstellar 
space or perturbed so that i t becomes a short-period comet 
ijWieeert fe Tremaind fl999l ). From the perspective of the 
Oort cloud, the comets are removed from the Oort cloud 
and from the loss co ne to the orbital parameter distribu- 
tion (|Fernandej|l98lf ). The loss cone (lc) is the population 
of comets that have orbits that will allow them to pene- 
trate the planeta ry system, making it possible to observe 
them (Hills 1981). Planetary perturbations move comets ef- 
ficiently from q values less than q\ c ~ 15 AU into either hy- 
perbolic orbits or into orbits that are more tightly bound to 
the solar system l|Hilldll98ll ). In the steady state situation, 
the new comets are distributed uniformly to the perihelion 
distances of q ^ gi c , for a > 30000 AU, while q > gi c comets 
come also from the inner region. For this reason, in all of 
our simulations we have assumed that initial perihelion dis- 
tances are distributed uniformly outside the loss cone. The 
inclination distribution of the outer Oort cloud of comets is 
isotropic (uniform in cosi). All the other angular elements 
are uniformly distributed between [0 — 27r]. 
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Figure 5. Variation of the parameters G\ (top panel), G2 (centre 
panel), and G3 (bottom panel) along the Solar orbit in the 'dy- 
namic' case (i.e. including full vertical and radial motions) over 1 
Gyr. 



3.1 Structure of the Oort cloud 

An important issue in studying the structure of the Oort 
cloud, is what energy distribution to adopt for the comets. 
We assume that the density of comets between 3000 AU and 
50000 AU in the Oort cloud is proportional to i/ r s - 5±0 - 5 } 
so that the n umber of comets N is dN oc l/r a dr, where 
a = 1.5 ± 0.5 l|Duncan et al.| [l987). The conclusions of this 
paper are not particularly sensitive to the adopted value of 
a. The existence of an inner Oort cloud has been speculated 
upon in many studies although there is no direct evidence for 
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Figure 6. Galactoccntric radius R, vertical height z and G3 for 
the 'dynamic' case of the Solar motion over 1 Gyr in the Galactic 
adopted model. Upper panel: the radial motion. Middle panel: 
vertical motion. Lower panel: G3. 



it (lHillslll98ll ). An inner Oort cloud is the extension of the 
Kuiper belt, filling the gap between the Kuiper belt and the 
outer Oort cloud. Semi-major axes i n the inner Oort cloud 
are typically between 50-15000 AU l|Leto et al.ll2008l ). The 
steady state flux from the inner Oort cloud cannot be uni- 
formly distributed in perihelion distance, since the planetary 
perturbations move comets efficiently from q values less than 
qi c ~ 15 AU to either hyperbolic orbits or i nto orbits that 
are more tightly bound to the solar system (Hills 1981). In 
the steady state situation, the new comets come uniformly 
to perihelion distances q ^ qi c when a > 30000 AU, while 
the q > qic comets also come from the inner region. 

3.2 Simulation parameters of the comets 

Due to this complicated picture, we study different semi- 
major axes in separate simulations, and evaluate the effi- 
ciency of tidal injection in each simulation. We chose initial 
sample conditions for the Oort cloud comets, setting the 
semi-major axis to be 10000, 20000, 30000, 40000, 50000, 
and 60000 AU. We choose the comet's eccentricity (e) ran- 
domly, so that the resulting values of q would be uniformly 
distributed from 35 to the value of the semi-major axis (a). 
All the other parameters, cos(i), Q, u, and the initial ec- 
centric anomaly were all chosen randomly in appropriate 
intervals. The number of comets in each simulation is 10 6 at 
each of the sampled semi-major axes. For the computation 
of the numbers of comets reaching the inner Solar System, 
the results from each semi-major axis are normalised to the 
adopted number density law. 

3.3 Simulations of the Galactic tide 

To study the effect of the Galactic potential on cometary 
motion, we chose two hypothetical solar systems: a 'con- 
stant' background density, where the Sun would be on a 
pure circular orbit with no vertical motion, and a realistic 
'dynamic' Solar orbit. In all systems, we analyse the flux 
of Oort cloud comets into the Solar System. This means 
that we consider a comet to have been detected in the in- 
ner Solar System when its q is within 30 AU, and it has 
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Table 3. Relative fraction of comets entering the inner Solar 
System, from each of the simulated semi-major axis values, for 
the 'constant' and 'dynamic' cases. 



Semi-major axis (AU) 


constant 


dynamic 


10000 


0.0729 


0.0732 


20000 


0.1803 


0.1804 


30000 


0.2405 


0.2405 


40000 


0.2381 


0.2369 


50000 


0.1504 


0.1507 


60000 


0.1179 


0.1184 



a heliocentric radius of less than 1000 AU. The last crite- 
rion is important, as the osculating elements of comets can 
evolve to have q ^ 30 AU, far away from the Sun, and evolve 
to more than 30 AU, without ever entering the inner Solar 
System. Since the minimum of the osculating q is at the 
aphelion, the criterion q ^ 30 AU is approximate. For this 
reason, we also check if the comet is actually approaching 
the aphelion, by choosing r < 1000 AU. This requires use of 
the mean anomaly, which is approximately calculated from 
the derived time co-ordinate. Comets which have been de- 
tected are removed from the simulation. We do not replace 
them with new comets. 



4 RESULTS 

Our aim is to determine the effect of the Galactic tide on 
Oort cloud comets by separately adopting a constant local 
density and a varying local density. Examining the flux of 
comets into the inner Solar System (Table [3]), we find that 
there is no significant difference between the two models. 
Most of the detected comets come from the middle (30—40 
kAU) range of semi-major axes. Likewise, the resulting dis- 
tribution of orbital elements, for comets coming into the 
Solar System, is the same in both the constant and dynamic 
cases. Fig. [7] shows example distributions for the simulation 
cases where the semi-major axis is 30000 AU. The gap at 
uj = and 180 is caused by the sin2aj 9 term going to in 
Equation [T] Similarly, the factor 1 — (J 2 /J 2 ) goes to zero 
at i=0 or 180, since J z = Jcosi. The g-distribution of in- 
coming comets shows no significant difference between the 
two cases, as seen in Fig. [8] The high peak at q = 30 AU, in 
Fig- [His caused by the slow evolution of q, at a = 30000 AU. 
Finally, the total flux of comets entering the inner-Solar sys- 
tem is not significantly different between the two cases. For 
the 'constant' case, from a source pool of 1.8 million comets 
in the simulated Oort cloud, we find that approximately 120 
comets per Myr reach the inner Solar System (which here 
means comets with q < 30 AU). Ass uming that there are 
abou t 10 12 comets in the Oort cloud (|Wiegert fc Tremaind 
1999), this corresponds to about 70 comets per year with 
q < 30 AU. This is a bit higher than the observed num- 
ber of new comets with q < 30 AU, which is about 20 per 
year, assuming a cometary flux of 0.65 ± 0.18 yr _1 AU - 
l|Fernandej|2010l ). 




i (degrees) 



Q (degrees) 





w (degrees) 

Figure 7. Distribution of the orbital elements i, Q, and w for 
a = 30000. The solid column represents constant case, and the 
shaded column the dynamic case. There is no significant difference 
between the two distributions. 




Figure 8. Distribution of the perihelion distance q for the in- 
coming comets, for a = 30000. The solid column represents the 
constant case, and the shaded column the dynamic case. There is 
no significant difference between the two distributions. 
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4.1 Time evolution of the cometary flux 

There is a clear difference between the two models when we 
look at the temporal evolution of the cometary flux. The 
top panel of Fig. shows, for the 'dynamic' case (and after 
a relaxation time of about 100 Myr) that the mean cometary 
flux (shown as a 10 Myr moving average) is well correlated 
with the changes in Gz (dashed line). The resulting fluxes 
from the simulation have been weighted according to their 
relative number-density (a -1,5 ), the resulting total amount 
of comets in the simulation, by using this weighting system 
is 1.8 xlO 6 . From Section r2.2l there are two main causes for 
the evolution of the flux, the major being the radial motion, 
the minor being the vertical motion. The top panel of Fig [9] 
shows the major trend clearly following the radial motion. 
The vertical motion is also followed, as is clear in the 10 
Myr moving average (which smooths out some of the Pois- 
son noise in the individual 1 Myr samples). We also ran a 
separate simulation with a circular orbit, and with the verti- 
cal component intact, the resulting fluxes followed perfectly 
the vertical motion, as was expected due to the changes in 
Gz being solely contributed by the vertical component. In 
the case where the values for Gi have been kept constant 
(i.e. corresponding to a completely circular orbit with no 
vertical motion), there is no evidence of any evolution, as 
seen in the bottom panel of Fig. [9] 

Simple Fourier-analysis of the dynamic flux (in the 
top panel of Fig. [9| finds two distinct periods in the 
cometary flux. The strongest signal is produced by a pe- 
riod of 143—167 Myr, and an equal signal from a period of 
41—45 Myr. The former corresponds to the radial period 
(152 Myr) and the latter to the half-period of the vertical 
period (43 Myr). The vertical signal is found to lie in the 
range 41 —45 Myr, and is qua si-periodic, as has been seen 
earlier bv lMatese et all (120011 ). 



4.2 Comparison with earlier results 

Matese et all (| 19951 ) derived various periods for the verti- 
cal Solar oscillation, depending on the adopted model of 
the disc. For example, their 'No-Dark-Disk Model' has a 
crossing-period very close to ours (43 Myr), while their 
'Best-fit Model' produces a much lower period of 33 Myr. 
They assumed a VK-velocity of 7.5 km s _1 , compared to 
ours of 7.25 km s _1 , and does not cause much qualitative 
difference in the vertical period of the orbit of the Sun. The 
lower period comes from assuming a considerably higher 
mid- plane density (p ~ 0.13 Mq/pc 3 ) than ours. 

iMatese et all l|200ll) found that the radial period of So- 
lar motion modulates the vertical period. However, they 
used very high local densities of matter in the disc, so 
that the vertical oscillations dominated the radial oscil- 
lations. This meant that the flux of comets into the in- 
ner Solar System with each passage through the disc was 
greatly amplified compared to our simulations. We find that 
the com etary flux var i es wit h an amplitude of about 20%, 
whereas IMatese et all (|200ll ) find flux variations of about 
a factor of two. The much smaller amplitude in the signal 
which we advocate, even taking into account the radial oscil- 
lations, would make finding a period in the scant cratering 
reco rd very difficult . 

iFouchard et al.l (|2006l ) calculated that a local mass den- 
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Figure 9. The time evolution of the relative flux of comets in the 
dynamic case (top panel), and the constant case (bottom panel). 
The solid line represents the relative flux of comets in 1 Myr bins, 
the white line a 10 Myr moving average. The dashed line in the 
top panel represents the evolution of the G3-parameter. 



sity of p = 0.1 Mq/pc 3 corresponds to a cometary flux of 
around 10 4 for their first 500 Myr interval for a source pop- 
ulation of 10 6 comets, where observed comets have q ^ 15 
AU. Assuming that there are 10 12 comets in the Oort cloud. 
This produces a flux of ~20 comets/yr, which is a factor of 
two less than our estimate. 

The Galactic tide case in I Packman et al l (|200ct) gives 
a flux of 100 comets per 50 Myr, from a source popula- 
tion of 10 6 comets, where observed comets have evolved 
from q > 15AU to q < 5 AU. This corresponds to a flux 
of 2 comets/yr, again assuming an Oort cloud with 10 12 
comets. Using similar analysis we find that we get a flux of 
290 comets per 50 Myr, correspo nding to a factor of three 
larger than lRickman et all (1200a). 

A review by iBailer- Jonesl (|2009l ) correlates terrestrial 
events with cometary signals. One of the more interesting 
ones in the revi e w is the 140 ± 15 Myr period found by 
iRohde fc Muiierl (2005) in the number of known marine an- 
imal genera as a function of time . While this c ould corre- 
spond to the Sun's radial period, Bailer- Jones (2009) con- 
siders that it has not been significantly detected, the main 
problem being that the entire time-span of the data covers 
no more than three oscillation s. Many other perio ds, proxies, 
and studies are mentioned in IBailer- Jonesl (|2009l ). although 
the conclusion is that there is no proven impact on biodi- 
versity as a result of the orbital motion of the Sun. 



5 DISCUSSION AND CONCLUSIONS 

We have studied the long-term dynamics of Oort cloud 
comets under the influence of both the radial and the ver- 
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tical components of the Galactic tidal field. Other perturb- 
ing forces on the comets, such as passing stars or passage 
through spiral arms are ignored, since we aim to study the 
influence of just the axisymmetric Galactic tidal field on the 
cometary motion. 

We use an axisymmetric model of the Galaxy, and a re- 
cently revised value for the Solar motion bv lSchonrich et all 
(2010). This leads to vertical oscillations of the Sun with 
an amplitude of about 100 pc, and radial oscillations over 
about 1 kpc. The changing tidal forces on the Oort cloud 
are computed as the Sun orbits for 1 Gyr in this potential, 
and the flux of comets entering the inner Solar System is 
computed in simulations. 

As expected, the cometary flux is strongly coupled to 
the G3-parameter in the tidal forces, which is dominated by 
the local mass density seen along the Solar orbit. Both the 
radial and vertical motions of the Sun can be seen in the 
cometary flux, although the amplitude of the variations is 
small, implying that detecting such a signal from the small 
number of age-dated craters would be very difficult. This 
agrees with the recent review of the detectability of the Solar 
motion in terrestrial proxies ( Bailer- Jones 2009). 

As Gs is directly coupled to local density, it is easily af- 
fected by the non-axisymmetric components in the motion 
of the Sun in the Galaxy. This implies that spiral arms, a 
Galactic bar, giant molecular clouds, or any other intermit- 
tently encountered structure should have an effect on the 
flux. 
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